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Abstract. Strange quark matter in a color flavor locked (CFL) state can be the true ground state of hadronic 
matter for a much wider range of the parameters of the model (the gap of the QCD Cooper pairs A, the strange 
quark mass m s and the Bag Constant B) than the state without any pairing. We review the equation of state 
(EOS) of CFL strange matter and study the structure of stellar objects made up of this phase, highlighting the 
novel features of the latter. Although the effects of pairing on the equation of state are thought to be small, we 
find that CFL stars may be in fact much more compact than strange stars (SS). This feature may be relevant 
in view of some recent observations claiming the existence of exotic and/or deconfined phases in some nearby 
neutron stars (NS). 
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1. Introduction 

The study of a hypothetical stability of strange quark mat¬ 
ter (SQM), put forward in Witten’s 1119841 seminal paper 
and a few important precursors (Bodmer 119711 Terazawa 
H3731 Chin and Kerman ri979l) has entered its third decade. 
Nothing less than the nature of the true ground state of 
hadronic matter is being questioned and in fact, within 
simple models, the hypothesis of a stable form of cold 
catalyzed plasma was shown to be tenable. Following the 
pioneering works, a general calculation of strange matter 
by Far hi and Jaffe in the framework of the MIT 

Bag model of confinement (Degrand et al. 119751 Cleymans, 
Suhonen and Gavai )15Sfi ) identified the so-called ’’win¬ 
dows of stability”, or regions in the plane m s — B in¬ 
side which the stability of SQM can be realized. Other 
models of confinement have been worked out to find a 
quite ample range of conditions for SQM to be absolutely 
bound (Zhang and Ru-Ken Su l2()U2[ and references therein; 
Hanauske et al. EED- However, there is a consensus that 
the issue of the availability of a ~ 1 % binding energy dif¬ 
ference for SQM to be bound is ultimately an experimental 
matter. 

While the search for SQM in laboratory and astro- 
physical environments beyond their present limits contin¬ 
ues, important theoretical developments have taken place. 
The most sound is probably the revival of interest in pair¬ 
ing interactions in dense matter, a subject already ad¬ 
dressed in the early 1980s (Bailin and Love 1198411 and 
revived a few years ago with new calculations of the pair¬ 


ing energy and related physics. It is now generally agreed 
(Alford, Rajagopal and Wilczek fl999bl Rapp et al. 120001 
Rajagopal and Wilczeck 120001 ) that (at least for very high 
densities) the color-flavor locked (CFL) state, with equal 
numbers of u, d and s quarks is likely to be the ground 
state of strong interactions, even if the quark masses are 
unequal (Alford, Berges and Raiagopal ll999al Schafer and 
Wilczek ll99 9l The equal number of flavors is enforced by 
the symmetry of the state, and thus electrons are absent 
because the mixture is automatically neutral (Rajagopal 
and Wilczeck [20011 Steiner, Reddy and Prakash 2002). 

Long awaited by theoretical physicists, the high- 
performance of the space X-ray missions Chandra and 
XMM (Weisskopf l2002l Becker and Aschenbach 12002 ) en¬ 
abled unprecedented studies of imaging and spectra of 
selected neutron stars with the aim of determining the 
masses and radii, perhaps the most simple forms of (in¬ 
directly) investigating the nature of high density mat¬ 
ter. Adopting General Relativity as a framework, a com¬ 
parison of the static models generated by integration of 
the Tolman-Oppenheimer-Volkof equations with observed 
data is expected to give information about the equation 
of state, and possibly other effects such as rotation, mag¬ 
netic atmospheres and so on. A wonderful example of the 
pre-space determinations is the mass of the binary pul¬ 
sar PSR 1913+16, accurate to several decimal places (van 
Kerkwiik l2DUl ). Other methods based on combinations of 
spectroscopic and photometric techniques have been re¬ 
cently proposed. At least one X-ray source is significantly 
above the centroid of the binary distribution ~ 1.4 M©; 
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namely Vela X-l for which a value of 1.87 Iq^Mq has 
been obtained (van KerkwijkjJDjllJ. Recently, not only the 
masses have been determined, but also indications of the 
radii became available, suggesting a very compact struc¬ 
ture. For instance, claims of high compactness have been 
made from the analysis of the binary Her X-l (Li, Lai and 
Wane ll995t Dey et al. 119980 with M = 0.98±0.12M e and 
R = 6.7 ± 1.2 km, and of the isolated nearby RX J185635- 
3754 (Pons et al. 12002 1 with M « 0.9Mg and i? « 6 km. 
In both cases, the results have been revisited and chal¬ 
lenged by other groups (Reynolds et al.! 19971 Kaplan, van 
Kerkwijk and Anderson IM who in turn found figures 
around the expected for conventional neutron star mod¬ 
els. This stresses the cautionary remarks made by several 
researchers about the high-compactness objects and guar¬ 
antees further studies, already undertaken in most cases. 
Li et al. ( 4999b ) have also added the source 4U 1728-34 
to the candidate list, showing that conventional accretion 
models indicate a very compact source in the mass-radius 
plane. As always, the actual distance to the source is a 
matter of concern. Needless to say, these results have yet 
to be confirmed carefully. Nonetheless, it is worthwhile to 
consider the possibility that at least some compact stars 
are extremely compact, or in other words, that their radii 
are ~ 30 — 40% less than the ’’canonical” lOfcm favored 
by neutron matter models for M ~ 1 Mq. 

If actually present in these sources, compactness would 
be extremely difficult (perhaps impossible) to model using 
underlying equations of state based on ordinary hadrons 
alone, and a natural alternative would be to consider de- 
confined matter (or other exotic components, like kaon 
matter or hyperons). Stars suspected to be made of decon- 
fined matter also include the X-ray bursters GRO J1744- 
28 (Cheng et al. 11998 ) and SAX J1808.4-3658 (Li at al. 
Il999at ). This increasing availability of data on compact 
neutron star structure coupled with the new developments 
on the QCD ground sate leads to investigate the proper¬ 
ties of compact stars in the light of the CFL state 1 . 

2. Color flavor locked strange matter 


sities. To proceed we need the thermodynamical potential 
Qcfl to derive the relevant quantities. Fortunately, this 
potential can be found quite simply (Alford et al. 1200 1 i to 
the lowest order in the gap A 2 . The procedure is to start 
with the f 1 f ree of a fictional state of unpaired quark mat¬ 
ter in which all quarks which are ’’going to pair” have a 
common Fermi momentum i/, with v chosen to minimize 
flfree- The binding energy of the diquark condensate is 
included by subtracting the condensation term 3A 2 /i 2 /7r 2 
while the vacuum energy is introduced by means of the 
phenomenological bag constant B allowing the mixture to 
confine. In this work we assume that the CFL state has 
been reached by deconfined matter, although the precise 
way in which this is done (involving weak decays, etc) 
merits a closer look. The thermodynamic potential f Icfl 
of this model reads (Alford et al. i20t)ll ) 


VlcFL = Qfree -y A 2 p 2 + B 

7 T 


.ZA 4^2 

i=u,d,s 


^ - lm^) + lmt log 


- ? A V + R, 


(1) 


where 3/t = + Hd, + Ms, Mi being the chemical potential 

of the i-species. The common Fermi momentum v = (/z 2 — 
m 2 ) 1 / 2 is given by 


V — 2 — I /i H- 


m : 


1/2 


( 2 ) 


The thermodynamic quantities, such as the pressure P, the 
baryon number density n B , the particle number densities 
rii , and the energy density e can be easily derived at T = 0 
and read (Lugones and Horvath 


P = —fi CFLy 


( 3 ) 


TlB — flu — Tld — Tls 


(v 3 + 2A 2 M ) 


( 4 ) 


2.1. The equation of state 

The very complex structure of the QCD phase diagram 
emerging from detailed calculations by several groups is 
still being explored. Therefore, only schematic models are 
available to explore stellar structure questions. It is, how¬ 
ever, widely agreed that if the quark mass m s is small 
enough (possibly the actual case in nature), the CFL state 
would be the minimum energy configuration at high den¬ 

1 Alford and Reddy (2002) have also independently analyzed 
the role of CFL quark matter in compact stars in a very recent 
work. Their work is rather general and focused on stars with a 
normal matter envelope, while our work deals with absolutely 

stable CFL stars only. However, it is worth noting that Alford 

and Reddy (2002) also include some results on pure CFL stars 
in their Figs. 5 and 6, which are in quantitative agreement with 
ours (see Section 3). 


£ = ^2 ^cfl = 3 fJ,n B - P. (5) 

i 

In this approach we shall treat the values of B, m s and 
A (possibly as high as ~ 100 MeV) as free constant pa¬ 
rameters. The full dependence of flcFL on m s has been 
known for years and certainly complicates the evaluation 
of the equation of state, which must be treated numeri¬ 
cally. However, we have recently worked out an approxi¬ 
mation to the order to 2 which has the main advantage of 
keeping the equation of state very simple, yet useful for 
most calculations, while at the same time highlights the 
effect of each parameter of the model (see Lugones and 
Horvath 120021 for details). 

The stiffness of the EOS is relevant for stellar models 
since it has a direct impact on the compactness of the stel¬ 
lar configurations. As discussed by Lugones and Horvath 
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Fig. 1. The dashed lines show the equation of state for 
CFL strange matter for B=75 MeV fm -3 , m s = 150 MeV, 
and two different values of the gap A as indicated in the 
figure. The solid line corresponds to SQM without pairing. 
Note the change of stiffness according to the value of A, 
as discussed in the text. 



Fig. 2. The stability windows for CFL strange matter. 
If the strange quark mass m s and the bag constant B 
lie inside the bounded region the CFL state is absolutely 
stable. Each window corresponds to a given value of the 
gap A as indicated by the label. The vertical solid line 
is the limit imposed by requiring instability of two-flavor 
quark matter. Note the enlargement of the window with 
increasing A. 


II2002I1 . whenever A is higher than m s / 2, the EOS is stiffer 
than the unpaired SQM (that is, produces more pressure 
for a given energy density). Since the actual value of A is 
not well known, we may expect either case, a stiffer or a 
softer EOS (for a given B ). 

2.2. Stability of the CFL phase 

Self-bound matter, in which strong forces produce a zero- 
pressure point at finite density, is potentially interesting 
for astrophysics because it produces stellar sequences hav¬ 


ing M oc R 3 in the Newtonian limit. For a given EOS 
the energy per baryon e/riB of the deconfined phase (at 
P = 0 and T = 0) must be lower than the neutron mass 
m n if matter is to be absolutely stable. Another condition 
to be considered results from the empirically known sta¬ 
bility of normal nuclear matter against deconfinement at 
P = 0 and T = 0 (Farhi and .Taffe 11984f) . which means 
that the energy per baryon of deconfined matter (a pure 
gas of quarks u and d ) at zero pressure and temperature 
must be higher than m n (neglecting refinements of order 
~ 1 MeV to the latter statement). Work using the MIT- 
based EOS has shown that the latter condition imposes 
that the MIT bag constant B must be larger than 57 MeV 
fm -3 (Farhi and .lalfe '. H)S I •. 

Since all three flavors have the same common Fermi 
momentum, the absolute stability condition takes a very 
simple form 


e 


n B 


p =o 


3/r < m n = 939 MeV. 


( 6 ) 


Therefore, at the zero pressure point we have (without any 
approximation) 


B — Flfree , Mo) V 9 A Poi 
7r 


(7) 


where po = m n /3 = 313 MeV. A family of curves in 
the m s — B plane follow, and on each one the energy per 
baryon is exactly e/n B = m n for a given A (Lugones and 
Horvath 2002). 

Fig. 2 displays the stability window for the CFL phase 
(i.e. the region in the m s versus B plane where E/n B is 
lower than 939 MeV at zero pressure. Eq. 0 sets the right 
side boundary of the window while the left side boundary 
is imposed by the minimum value B = 57 MeV. The win¬ 
dow is greatly enlarged for increasing values of A with 
respect to the original SQM calculations which do not in¬ 
clude any pairing (see for example, Fig. 1 of Farhi and 
.Taffe llOffl] ). We have emphasized elsewhere that the pair¬ 
ing gap A may be important for a ” CFL strange matter” 
state (Lugones and Horvath 2002, see also Madsen ;200l! 
for a similar view). 


3. Structure of color flavor locked stars 

The study of the nrass-radius relation for compact stars is 
a widely known tool for testing the existence of different 
phases of matter inside NSs. It is also well-known that si¬ 
multaneous observational data on masses and radii can im¬ 
pose important constraints on the high density equations 
of state. Up to now, almost all of the measured masses of 
’’neutron stars” clustered within a narrow range around 
1.4 Mg. It has been conjectured that this mass scale may 
be due to the fact that neutron stars are formed in the 
gravitational collapse of supernovae and come just from 
the iron cores of the pre-supernovae. However, nothing 
fundamental precludes smaller neutron stars from exist¬ 
ing provided some mechanism capable of creating them 
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Fig. 3. The mass-radius relation for CFL strange stars for 
B =115 MeV fm -3 , m s = 150 MeV, A = 100 MeV on the 
left, and for B = 70 MeV fm -3 , m s = 150 MeV, A = 
100 MeV on the right. The crosses indicate the models 
of maximum radius for each set of parameters, while the 
asterisks indicate the models with the maximum mass (see 
Figs 4-7 for more details). 

operates. Because small masses are not expected in core 
collapse models, they hold the potential for discriminat¬ 
ing among possible equations of state. This is particularly 
significant given the recent claims of very low-mass/low- 
radii objects. As already mentioned, at least two sources, 
Her X-l and RX J1856-37, are candidates for high com¬ 
pactness. They have been claimed to have radii ~ 7 km 
and masses around 1M 0 (Li et al. 11995 : Pons et al. K(lf)2(l . 
Fragmentation of a larger star (perhaps due to high ro¬ 
tation) may be a possible way to produce sub-solar mass 
compact objects ('Nakamura l2()()2ll . Nevertheless, it is clear 
that models and data are still subject to detailed analysis. 

CFL strange matter offers the double bonus of enhanc¬ 
ing the odds for a self-bound state, and also of allow¬ 
ing the existence of stable stars made up entirely of this 
phase. To pinpoint the effects of such a composition we 
have solved the Tolman-Oppenheimer-Volkoff equations of 
stellar structure, assuming a parameter set inside the win¬ 
dows of absolute stability discussed above (Lugones and 
Horvath mm- Fig. 3 shows sequences of compact star 
models thus calculated for particular values of the param¬ 
eters A, m s and B (see caption). These values of the latter 
quantities have been selected to highlight a potentially im¬ 
portant feature of pairing in dense quark matter, namely 
the possibility of obtaining very compact stellar models. 
As expected, the mass-radius relation of stars made up of 
CFL strange matter presents the same qualitative shape 
as the curves found for strange stars . This is a direct con¬ 
sequence of the existence of a zero pressure point at finite 
density (of the order of the nuclear saturation density). 

To assess the robustness of the results, and to under¬ 
stand the outcome of masses and radii for a given set of 
parameters, we have explored the effect of each parameters 
of the EOS on the static structure of the stellar models. 


As shown in Figs. 4-7, and quite analogously to the case 
of strange stars, the effect of increasing B while holding A 
and m s fixed is to lower both the maximum radius R ma x 
and the maximum mass M max of the stellar configura¬ 
tions. The mass of the strange quark m s works in the same 
direction although the effect is much smaller. The effect of 
the pairing gap A is qualitatively different and much more 
interesting. An increase of A increases the value of both 
M max and R m ax- However, since the increase of A strongly 
enlarges the stability window, the combination of a large 
A with a large B allows the construction of very compact 
models made up of absolutely stable CFL strange matter. 
This is particularly illustrative when we set A = m s / 2. 
In this case the EOS adopts, to order /.t 2 , the same sim¬ 
ple form e = 3 P + 4 B as unpaired SQM. However, since 
the stability window is considerably enlarged, large val¬ 
ues of B are allowed and very compact stars are found. 
This shows how, although the effects in the EOS are only 
a few percent (O ( A/fj ,) 2 ), the enlargement of the stabil¬ 
ity window allows a very rich set of stellar configurations, 
ranging from very compact to quite extended objects. As 
is seen from the Figures, for the smallest values of B the 
models tend to be less compact and show radii up to 13 
km and masses as high as 2.4 M 0 . The most extended 
configurations are obtained for a combination of small B 
and high A. 2 It is not possible to select the most likely 
range of the parameters, which should be limited by ob¬ 
servational arguments much in the same way as done for 
strange matter models. 

4. Discussion 

We have addressed in this work stellar models constructed 
with the simplest version of a ”CFL strange matter” 
phase. The CFL phase at zero temperature has been mod¬ 
elled as an electrically neutral and colorless gas of quark 
Cooper pairs, in which quarks are paired in such a way 
that all the flavors have the same Fermi momentum and 
hence the same number density (Rajagopal and Wilczeck 
mm- Since the strange quark is massive, some energy 
must be spent in order to keep a common Fermi mo¬ 
mentum for all three flavors. However, more than that 
is gained from the energy gap of the pairing. The model 
allows CFL strange matter to be the true ground state of 
strong interactions for a wide range of the parameters B , 
m s and A. It is remarkable that the strange matter hy¬ 
pothesis may be quite boosted by pairing interactions. We 
have explored the mass-radius relation for all the values of 

2 It can be checked that when the two particular values of 

B used by Alford and Reddy (2002) are selected (B 1 / 4 = 165 
MeV and B 1/,4 = 185 MeV; corresponding to B = 97 MeV fm -3 

and B = 153 MeV fm~ 3 respectively) the values obtained for 
Mmax and Rmax from Figs. 4-7 are quantitatively in agree¬ 
ment (even if the latter value has to be extrapolated if needed 
in these Figs.). Moreover, a variety of CFL strange star mod¬ 
els can be explored by varying the parameters, although, for 
example, a large variation in m s does not result in substantial 
modifications of Mmax and R ma x- 
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Fig. 4. The maximum radius indicated with crosses in Fig. 
3 is shown here as a function of B, for m s = 150 MeV and 
different values of A. 



Fig. 6. The maximum mass indicated with asterisks in 
Fig. 3 is shown here as a function of B, for m s = 150 MeV 
and different values of A. 



Fig. 5. The same as the previous figure but for fixed A 
and different values of m s . 



Fig. 7. The same as the previous figure but for fixed A 
and different values of m s . 


the parameter space of the EOS that give absolutely sta¬ 
ble CFL strange matter. Very compact configurations are 
found that could help to explain the recently claimed com¬ 
pactness of a few neutron stars. Also, for some parameter 
choices the models are massive and extended. 

A feature that has been largely explored and debated 
in the context of strange star models is the existence of 
a normal matter crust (Zdunik [200211 . Bare quark sur¬ 
faces may alter drastically the ability of radiating pho¬ 
tons, and thus they may produce interesting signatures 
for their identification (Page and Usov 120(121) . Conversely, 
a normal matter crust (held in mechanical equilibrium by 
the electrostatic potential at the surface) may hide most 
of the features of exotic matter. By its very construction, 
it is clear that (in striking contrast with the well-studied 
SS) no crust could be present in the case of CFL strange 
stars. This is directly related to the absence of electrons in 
the mixture (Rajagopal and Wilczek 120011 . and thus to 
the absence of an electrostatic potential to prevent nor¬ 
mal matter from being converted to the stable CFL state. 


CFL strange stars must have bare surfaces within these 
models. This in turn means that the transport properties 
of the state, which are still being explored (see, for ex¬ 
ample Shovkovy and Ellis 1200'S : Jaikumar, Prakash and 
Schafer 2002; Reddy, Sadzikowski and Tachibana 2003), 
hold the clue to the identification of these compact stars. 
We expect the photon emission properties of CFL strange 
stars to resemble those of bare SS. Briefly, pairing effects 
should affect the plasma frequency oj p through the baryon 
number density as a small correction of order /rA 2 so that 
uj p will not be very different from the ~ 20 MeV expected 
for bare SS . The equilibrium photon radiation will show 
a very hard spectrum and a tiny luminosity, making CFL 
strange stars very difficult to detect. On the other hand, 
while the thermal emission of photons from the bare quark 
surface of a hot strange star (mainly by electron-positron 
pair production) has been shown to be much higher than 
the Eddington limit (Page and UsovfSOiEJ, this would not 
be the case for CFL SS since no electrons will be present 
at the surface. 
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As already mentioned, the fact of considering CFL 
SQM to be absolutely stable together with the absence of 
electrons in the mixture precludes the existence of a crust 
of normal matter, then CFL strange stars are bare by con¬ 
struction. In a more general approach, Alford and Reddy 
(2002b) found mass-radius relations with essentially the 
same shape as that for neutron stars made up of normal 
nuclear matter when considering a wider range of parame¬ 
ters in which ’’hybrid” stars are allowed. This corresponds 
to a large part of the parameter space and not only do they 
find a normal matter envelope for all these models, but 
also the composition of the less massive objects becomes 
pure nuclear matter, as expected. In addition, as pointed 
out in Section 3, a subset of the models of Alford and 
Reddy (2002) correspond to our ”CFL strange stars”, and 
the full range of these self-bound models can be obtained 
with the aid of Figs. 4-7. The self-bound CFL state may 
be compared to stable diquark states, which have been 
found to produce very similar stellar structural proper¬ 
ties than CFL strange matter within completely different 
models (Horvath IlDi).*)! Horvath and de Freitas Pacheco 
llQ98 Lugones and Horvath 12003 L We finally note that, 
due to the presence of m s and A in the EOS, a linear 
form P(p) is obtained in particular cases only, therefore 
the scaling behavior of the Tolman-Oppenheimer-Volkov 
equations that allows the construction of simple expres¬ 
sions for the maximum masses and radii is lost for these 
CFL stars. The relevance of this modelling to understand 
actual observed stars is not clear as yet, but certainly the 
latter will continue to advance to the point at which this 
and other questions will be answered with confidence. 
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